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We discuss a new phenomenon of turbulent thermal diffusion associated with turbulent transport 
of aerosols in the atmosphere and in laboratory experiments. The essence of this phenomenon 
is the appearance of a nondiffusive mean flux of particles in the direction of the mean heat flux, 
which results in the formation of large-scale inhomogeneities in the spatial distribution of aerosols 
that accumulate in regions of minimum mean temperature of the surrounding fluid. This effect of 
turbulent thermal diffusion was detected experimentally. In experiments turbulence was generated 
by two oscillating grids in two directions of the imposed vertical mean temperature gradient. We 
used Particle Image Velocimetry to determine the turbulent velocity fleld, and an Image Processing 
Technique based on an analysis of the intensity of Mie scattering to determine the spatial distribution 
of aerosols. Analysis of the intensity of laser light Mie scattering by aerosols showed that aerosols 
accumulate in the vicinity of the minimum mean temperature due to the effect of turbulent thermal 
diffusion. 



I. INTRODUCTION 

Aerosols are a universal feature of the Earth's atmo- 
sphere. They can significantly affect the heat balance 
and dynamics of the atmosphere, climate, atmospheric 
chemistry, radiative transport and precipitation forma- 
tion (see, e.g., Twomey, 1977; Seinfeld, 1986; Flagan and 
Seinfeld, 1988; Pruppacher and Klett 1997; Lohmann and 
Lesins, 2002; Kaufman et al. 2002; Anderson et al. 2003; 
and references therein) . Formation of aerosol clouds is of 
fundamental significance in many areas of environmental 
sciences, physics of the atmosphere and meteorology (see, 
e.g., Seinfeld, 1986; Flagan and Seinfeld, 1988; Paluch 
and Baumgardner, 1989; Shaw et al. 1998; Shaw 2003; 
and references therein). It is well known that turbulence 
causes decay of inhomogeneities in spatial distribution of 
aerosols due to turbulent diffusion (see, e.g., Csanady, 
1980; McComb 1990; Stock 1996), whereas the opposite 
effect, the preferential concentration of aerosols in atmo- 
spheric turbulence is not well understood. 

Elperin et al. (1996; 1997; 1998; 2000b; 2001) recently 
found that in a turbulent fluid flow with a nonzero mean 
temperature gradient an additional mean flux of aerosols 
appears in the direction opposite to that of the mean tem- 
perature gradient, which is known as the phenomenon of 
turbulent thermal diffusion. This effect results in the 
formation of large-scale inhomogeneities in the spatial 
distribution of the aerosol particles. 

The phenomenon of turbulent thermal diffusion is im- 
portant for understanding atmospheric phenomena (e.g., 
atmospheric aerosols, smog formation, etc.). There ex- 
ists a correlation between the appearance of tempera- 



ture inversions and the aerosol layers (pollutants) in the 
vicinity of the temperature inversions (see, e.g., Seinfeld, 
1986; Flagan and Seinfeld, 1988). Moreover, turbulent 
thermal diffusion can cause the formation of large-scale 
aerosol layers in the vicinity of temperature inversions in 
atmospheric turbulence (Elperin et al. 2000a; 2000c). 

The main goal of this paper is to describe the experi- 
mental detection of a new phenomenon of turbulent ther- 
mal diffusion. To accomplish this, we constructed an 
oscillating grids turbulence generator (for details see Ei- 
delman et al. 2002). Recent studies by De Silva and 
Fernando (1994); Srdic et al. (1996), Shy et al. (1997) 
have demonstrated the feasibility of generating nearly 
isotropic turbulence by two oscillating grids. Turbulent 
diffusion in oscillating grids turbulence was investigated 
by Ott and Mann (2000). In order to study turbulent 
thermal diffusion we used Particle Image Velocimetry to 
characterize a turbulent velocity field, and an Image Pro- 
cessing Technique based on Mie scattering to determine 
the spatial distribution of the particles. Our experiments 
were performed in two directions of the vertical mean 
temperature gradient: an upward mean temperature gra- 
dient (formed by a cold bottom and a hot top wall of the 
chamber) and a downward mean temperature gradient. 
We found that in a flow with a downward mean temper- 
ature gradient, particles accumulate in the vicinity of the 
top wall of the chamber (whereby the mean fluid temper- 
ature is minimal), and in a flow with an upward mean 
temperature gradient particles accumulate in the vicin- 
ity of the bottom wall of the chamber due to the effect 
of turbulent thermal diffusion. 

The paper is organized as follows. Section II discusses 
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the physics of turbulent thermal diffusion, and Section III 
describes the experimental set-up for a laboratory study 
of this effect. The experimental results are presented 
in Section IV, and a detailed analysis of experimental 
detection of turbulent thermal diffusion is performed in 
Section V. Finally, conclusions arc drawn in Section VI. 



II. TURBULENT THERMAL DIFFUSION 

Evolution of the number density n{t, r) of small parti- 
cles in a turbulent flow is determined by 
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where J„ is the molecular flux of the particles and Vp 
is the velocity of the particles in the turbulent fluid ve- 
locity field. Averaging Eq. ^ over the statistics of the 
turbulent velocity field yields the following equation for 
the mean number density of particles N = (n) : 
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where = (r/3)(u'^) is the turbulent diffusion coeffi- 
cient, r is the momentum relaxation time of the turbu- 
lent velocity field, Vp — Vp -I- u, Vp — (vp) is the mean 
particle velocity, and the effective velocity is 



(ru(V • u)) . 



(4) 



Here J^,, = ~D{'VN + ktVT/T) is the mean molecular 
flux of particles, D is the coefficient of molecular diffu- 
sion, kt is the thermal diffusion ratio, and T = (T) is 
the mean fluid temperature. Equations ^ and were 
previously derived by different methods (see Elperin et 
al. 1996; 1997; 1998; 2000b; 2001; Pandya and Mashayek 
2002). 

For nonincrtial particles advected by a turbulent fluid 
flow, particle velocity Vp coincides with fluid velocity v, 
and V • V « -(v • V)p/p « (v • V)T/T, where p and 
T are the density and temperature of the fluid. Thus, 
the effective velocity (0} for nonincrtial particles can be 
given by Veff = —D^{WT)/T, which takes into account 
the equation of state for an ideal gas but neglects small 
gradients of the mean fluid pressure. 

For inertial particles, velocity Vp depends on the veloc- 
ity of the surrounding fluid v. Velocity Vp can be deter- 
mined by the equation of motion for a particle. Solving 
the equation of motion for a small solid particle with 
Pp 3> p yields: Vp = v — Tpdv/dt -\- 0(rp) (see Maxey 
1987), where Tp is the Stokes time and pp is the material 
density of the particles. In that case 



AP 

V • Vp = V • V Tp— + 0(t2) , 



(5) 



and the effective velocity of the inertial particles can be 
given by Vos = -D^{1 + k)(VT)/T, where P is the 



fluid pressure. Coefficient k depends on particle inertia 
{nip/m^), the parameters of turbulence (Reynolds num- 
ber) and the mean fluid temperature (Elperin et al. 1996; 
1997; 1998; 2000b; 2001). Here rrip is the particle mass 
and is the mass of molecules of the surrounding fluid. 
The turbulent flux of particles Q can be rewritten as 
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where = (1 -I- k)N. The flrst term in the RHS of 
Eq. © describes turbulent thermal diffusion, while the 
second term in the turbulent flux of particles ^ describes 
turbulent diffusion. Parameter fc^ can be interpreted as 
the turbulent thermal diffusion ratio, and D^k^ is the co- 
efficient of turbulent thermal diffusion. Turbulent ther- 
mal diffusion causes the formation of a large-scale pat- 
tern wherein the initial spatial distribution of particles 
in a turbulent fluid flow evolves into a large-scale inho- 
mogeneous distribution, i.e., particles accumulate in the 
vicinity of the minimum mean temperature of the sur- 
rounding fluid (Elperin et al. 1996; 1997; 1998; 2000b; 
2001). 

The mechanism responsible for the occurrence of tur- 
bulent thermal diffusion for particles with pp ^ p can be 
described as follows. Inertia causes particles inside the 
turbulent eddies to drift out to the boundary regions be- 
tween the eddies (i.e., regions with low vorticity or high 
strain rate and maximum of fluid pressure). Thus, parti- 
cles accumulate in regions with maximum pressure of the 
turbulent fluid. For simplicity, let us consider a pure in- 
ertial effect, i.e., we assume that V • v = 0. This inertial 
effect results in V • Vp esc TpAP 0. On the other hand, 
Eq. for large Peclet numbers yields V • Vp cx —dn/dt. 
The latter formula implies that dn/dt cx —TpAP, i.e., 
inertial particles accumulate [dn/dt > 0) in regions 
with maximum pressure of the turbulent fluid (where 
AP < 0). Similarly, there is an outflow of particles 
from regions with minimum pressure of fluid. In ho- 
mogeneous and isotropic turbulence without large-scale 
external gradients of temperature, a drift from regions 
with increased or decreased concentration of particles by 
a turbulent flow of fluid is equiprobable in all directions, 
and pressure and temperature of the surrounding fluid 
do not correlate with the turbulent velocity field. Thus 
only turbulent diffusion of particles takes place. 

In a turbulent fluid flow with a mean temperature gra- 
dient, the mean heat flux (uO) is not zero, i.e., the fluc- 
tuations of fluid temperature Q = T — T and the velocity 
of the fluid correlate. Fluctuations of temperature cause 
fluctuations of fluid pressure. These fluctuations result in 
fluctuations of the number density of particles. Indeed, 
an increase in pressure of the surrounding fluid is accom- 
panied by an accumulation of particles. Therefore, the 
direction of the mean flux of particles coincides with that 
of the heat flux, (vpn) c>c (u8) cx — VT, i.e., the mean 
flux of particles is directed to the area with minimum 
mean temperature, and the particles accumulate in that 
region. 
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FIG. 1: The experimental set-up. 




FIG. 2: The test section of the oscillating grids turbulence 
generator. 



III. EXPERIMENTAL SET-UP 

In this Section we investigate experimentally the effect 
of turbulent thermal diffusion. The experiments were 
conducted in an oscillating grids turbulence generator 
in air flow (see Fig. 1). The test section consisted of 
a rectangular chamber of dimensions 29 x 29 x 58 cm 
(see Fig. 2). Pairs of vertically oriented grids with bars 
arranged in a square array were attached to the right 
and left horizontal rods. Both grids were driven indepen- 
dently with speed-controlled motors. The grids were po- 
sitioned at a distance of two-grid meshes from the cham- 
ber walls parallel to them. A two-grid system can oscil- 
late at a controllable frequency up to 20 Hz. The grid 
stroke was adjusted within a range of 1 to 10 cm. 

A vertical mean temperature gradient in the turbulent 



flow was formed by attaching two aluminium heat ex- 
changers to the bottom and top walls of the test section 
(see Fig. 2). The experiments were performed in two 
directions of the mean temperature gradient: an upward 
mean temperature gradient (a cold bottom and a hot top 
wall of the chamber) and a downward mean temperature 
gradient (a heated bottom and a cold top wall of the 
chamber). In order to improve the heat transfer in the 
boundary layers at the walls we used a heat exchanger 
with rectangular pins 3 x 3 x 15 mm. This allowed us 
to support a mean temperature gradient in the core of 
the flow up to 200 K/m in the upward direction and up 
to 110 K/m in the downward direction at a mean tem- 
perature of about 300 K. The temperature was measured 
with a high-frequency response thermocouple. 

The velocity field was measured using a Particle Im- 
age Velocimetry (PIV), see Raffel et al. (1998). A digi- 
tal PIV system with LaVision Flow Master III was used. 
A double-pulsed light sheet was provided by a Nd-YAG 
laser source (Continuum Surelitc 2 x 170 mJ). The light 
sheet optics includes spherical and cylindrical Galilei tele- 
scopes with tuneable divergence and adjustable focus 
length. We used a progressive-scan 12 bit digital CCD 
camera (pixel size 6.7 /xm x6.7/im each) with a dual- 
frame-technique for cross-correlation processing of cap- 
tured images. A programmable Timing Unit (PC inter- 
face card) generated sequences of pulses to control the 
laser, camera and data acquisition rate. The software 
package DaVis 6 was applied to control all hardware com- 
ponents and for 32 bit image acquisition and visualiza- 
tion. This software package contains PIV software for 
calculating the flow flelds using cross-correlation analy- 
sis. Velocity maps and their characteristics, e.g., statis- 
tics and PDF, were analyzed with this package plus an 
additional developed software package. An incense smoke 
with sub-micron particles (with pp/ p ~ 10^) as a tracer 
was used for the PIV measurements. Smoke was pro- 
duced by high temperature sublimation of solid incense 
particles. Analysis of smoke particles using a microscope 
(Nikon, Epiphot with an amplification of 560) and a PM- 
300 portable laser particulate analyzer showed that these 
particles have a spherical shape and that their mean di- 
ameter is 0.7/im. 

Mean and r.m.s. velocities, two-point correlation func- 
tions and an integral scale of turbulence from the mea- 
sured velocity fields were determined. A series of 100 
pairs of images acquired with a frequency of 4 Hz were 
stored for calculating the velocity maps and for ensemble 
and spatial averaging of turbulence characteristics. The 
center of the measurement region coincides with the cen- 
ter of the chamber. We measured the velocity for flow 
areas from 60 x 60 mm^ up to 212 x 212 mm^ with a 
spatial resolution of 1024 x 1024 pixels. This size of the 
probed area corresponds to a spatial resolution from 58 
/Ltm / pixel up to 207 /xm / pixel. These regions were an- 
alyzed with interrogation windows of 32 x 32 and 16 x 16 
pixels. A velocity vector was determined in every in- 
terrogation window, allowing us to construct a velocity 
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map comprising 32 x 32 or 64 x 64 vectors. The veloc- 
ity maps were determined in two planes. In the one-grid 
experiments the plane was parallel to the grid, while in 
the two-grid experiments the plane was normal to the 
grids. The mean and r.m.s. velocities for each point of 
the velocity map (1024 points) were determined by av- 
eraging over 100 independent maps, and then over 1024 
points. The two-point correlation functions of the ve- 
locity field were determined for each point of the cen- 
tral part of the velocity map (16 x 16 vectors) by av- 
eraging over 100 independent velocity maps, and then 
over 256 points. An integral scale L of turbulence was 
determined from the two-point correlation functions of 
the velocity field. These measurements were repeated for 
different distances from the grid, and for various tem- 
perature gradients, Reynolds numbers and particle mass 
loadings. The PIV measurements performed in the oscil- 
lating grids turbulence generator confirmed earlier results 
(Thompson and Turner 1975; Hopfinger and Toly 1976; 
Kit ct al. 1997) for the dependence of varies character- 
istics of the turbulent velocity field on the parameters of 
the oscillating grids turbulence generator. 

To determine the spatial distribution of the particles, 
we used an Image Processing Technique based on Mie 
scattering. The advantages of this method have been 
demonstrated by Guibert et al. (2001). The light ra- 
diation energy flux scattered by small particles is Eg cx 
E '^{-Kdp/X; a; n), where E oc 7rdp/4 is the energy flux in- 
cident at the particle, dp is the particle diameter, A is the 
wavelength, a is the index of refraction and is the scat- 
tering function. In the general case, is given by Mie 
equations. For wavelengths A, which are smaller than the 
particle size, ^ tends to be independent of dp and A. The 
scattered light energy flux incident on the CCD camera 
probe is proportional to the particle number density n, 
i.e., Es oc En{ndpl4). 

Mie scattering does not change from temperature ef- 
fects since it depends on the permittivity of particles, the 
particle size and the laser light wavelength. The effect 
of the temperature on these characteristics is negligibly 
small. Note that in each experiment before taking the 
measurements, we let the system run for some time (up 
to 30 minutes after smoke injection into the flow with 
a steady mean temperature profile) in order to attain a 
stationary state. 

We found that the probability density function of the 
particle size measured with the PM300 particulate ana- 
lyzer was to the most extent independent of the location 
in the flow for incense particle size of 0.5 — 1 /xm. Note 
that since the number density of the particles is small 
(about 1 mm apart), it can be assumed that a change 
in particle number density will not affect their size dis- 
tribution. Therefore, the ratio of the scattered radiation 
fluxes at two points and at the image measured with the 
CCD camera remains equal to the ratio of the particle 
mean number densities at these two locations. 



IV. EXPERIMENTAL DETECTION OF 
TURBULENT THERMAL DIFFUSION 

The turbulent flow parameters in the oscillating grids 
turbulence generator are: r.m.s. velocity \J (u^) =4—14 
cm/s depending on the frequency of the grid oscillations, 
integral scale of turbulence L = 1.6 — 2.3 cm, and the 
Kolmogorov length scale 77 = 380 — 600 ^m. Other pa- 
rameters are given in Section 5 (Tables 1-2). These flows 
involve a wide range of scales of turbulent motions. In- 
terestingly, a flow with a wide range of spatial scales is 
already formed at comparatively low frequencies of grid 
oscillations. We found a weak mean flow in the form 
of two large toroidal structures parallel and adjacent to 
the grids. The interaction of these structures results in 
a symmetric mean flow that is sensitive to the param- 
eters of the grid adjustments. We particularly studied 
the parameters that affect a mean flow such as the grid 
distance to the walls of the chamber and partitions, and 
the angles between the grid planes and the axes of their 
oscillations. Varying these parameters allowed us to ex- 
pand the central region with a homogeneous turbulence. 
We found that the measured r.m.s. velocity was several 
times higher than the mean velocity in the core of the 
flow. 

The temperature measurements were taken in the os- 
cillating grids turbulence generator. The temperature 
difference between the heat exchangers varied within a 
range of 25 to 50 K. The mean temperature vertical pro- 
files at a frequency of grids oscillations / = 10.5 Hz in 
turbulent flows with downward and upward mean tem- 
perature gradients are shown in Fig. 3. Here Z is a di- 
mensionless vertical coordinate measured in units of the 
height of the chamber, and Z — Q dX the bottom of the 
chamber. Hereafter, we use the following system of coor- 
dinates: Z is the vertical axis, the F-axis is perpendicu- 
lar to the grids and the XZ-plane is parallel to the grids 
plane. 

Measurements performed using different concentra- 
tions of incense smoke showed that in an isothermal flow 
the distribution of the average scattered light intensity 
over a vertical coordinate is independent of the mean par- 
ticle number density. In order to characterize the spatial 
distribution of particle number density N oc E^ jE in 
a non-isothermal flow, the distribution of the scattered 
light intensity E for the isothermal case was used to nor- 
malize the scattered light intensity i?^ obtained in a non- 
isothermal flow under the same conditions. The scattered 
light intensities E^^ and E in each experiment were nor- 
malized by corresponding scattered light intensities av- 
eraged over the vertical coordinate. The ratios E'^ /E of 
the normalized average distributions of the intensity of 
scattered light as a function of the normalized vertical 
coordinate Z in turbulent flows with a downward mean 
temperature gradient and with an upward mean temper- 
ature gradient are shown in Fig. 4. The figure demon- 
strates that particles are redistributed in a turbiflent flow 
with a mean temperature gradient. Particles accumulate 
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FIG. 3: Vertical temperature profiles at frequency of grid 
oscillations / = 10.5 Hz in turbulent flows with an upward 
moan temperature gradient (filled squares) and with a down- 
ward mean temperature gradient (unfilled triangles). Here 
Z is a dimensionless vertical coordinate measured in units of 
the height of the chamber, and Z = at the bottom of the 
chamber. 



in regions of minimum mean temperature (in the lower 
part of the chamber in turbulent flows with an upward 
mean temperature gradient), and there is an outflow of 
particles from the upper part of the chamber where the 
mean temperature is larger. On the other hand, in a 
flow with a downward mean temperature gradient (a hot 
bottom and a cold top wall of the test section) particles 
accumulate in the vicinity of the top wall of the chamber, 
i.e., in the vicinity where the mean fluid temperature is 
minimum. 

To determine the turbulent thermal diffusion ratio, we 

plotted the normalized particle number density = 
N/No versus the normalized temperature difference Tr = 
(T — To) /To, where Tq is the reference temperature and 
No = N{f = To). The function Nr{Tr) is shown in Fig. 5 
for a turbulent flow with an upward mean temperature 
gradient and in Fig. 6 for a turbulent flow with a down- 
ward mean temperature gradient. 

We probed the central 20 x 20 cm region in the chamber 
by determining the mean intensity of scattered light in 

32 X 16 interrogation windows with a size of 32 x 64 pix- 
els. The vertical distribution of the intensity of the scat- 



FIG. 4: Ratios /E of normalized average distributions of 
the intensity of scattered light versus the normalized vertical 
coordinate Z in turbulent fiows with an upward mean tem- 
perature gradient (filled squares) and with a downward mean 
temperature gradient (unfilled triangles) . Here Z = is at the 
bottom of the chamber. The frequency of grids oscillations is 
/ = 10.5 Hz. 

tered light was determined in 16 vertical strips, which are 
composed of 32 interrogation windows. Variations of the 
obtained vertical distributions between these strips were 
very small. Therefore in Figs. 4-6 we used spatial aver- 
age over strips and the ensemble average over 100 images 
of the vertical distributions of the intensity of scattered 
light. Figures 5-6 were plotted using the mean temper- 
ature vertical profiles shown in Fig. 3. The normalized 
local mean temperatures [the relative temperature differ- 
ences {T — To) /To] in Figs. 5-6 correspond to the different 
points inside the probed region. In particular, in Fig. 5 
the location of the point with reference temperature Tq 
is Z = 0.135 (the lowest point of the probed region with 
a maximum N in turbulent flows with an upward mean 
temperature gradient), and the point with a maximum 
normalized temperature difference is the highest point of 
the probed region at Z = 0.865. The size of the probed 
region did not affect our results. 

V. DISCUSSION 

The vertical profiles of the mean number density of 
particles obtained in our experiments can be explained 
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FIG. 5: Normalized particle number density Nr = N/Nq 
versus normalized temperature difference Tr = {T — Tq) /Tq in 
a turbulent flow with an upward mean temperature gradient. 
The frequency of grid oscillations is / = 10.5 Hz. 



FIG. 6: Normalized particle number density Nr = N/Nq 
versus normalized temperature difference Tr = {T — Tq) /Tq in 
a turbulent flow with a downward mean temperature gradient. 
The frequency of grid oscillations is / = 10.4 Hz. 



by considering Eqs. and If one does not take into 
account the term TVVofr in Eq. jSJ for the mean number 
density of particles, the equation reads dN/dt = D^AN, 
where we neglected a small mean fluid velocity Vp and a 
small molecular mean flux of particles (that corresponds 
to the conditions of the experiment). The steady-state 
solution of this equation is A'^ = const. However, our 
measurements demonstrate that the solution N = const 
is valid only for an isothermal turbulent flow, i.e., when 
T = const. By taking into account the effect of turbu- 
lent thermal diffusion, i.e., the term NV^g in Eq. the 
steady-state solution of this equation for non-inertial par- 
ticles is: VN/N ~ -Vf/f. The latter equation yields 
N/No = 1 - (T - To) /To, where T-Tq-^ To- In the ex- 
periments using particles of the size 0.5 — 2 /im, wc found 
that 



where the coefficient a — 1.33 — 1.79 (depending on 
the frequency of the grid oscillations) in a turbulent 
flow with an upward mean temperature gradient, and 
a = 1.29—1.87 in a turbulent flow with a downward mean 
temperature gradient (see Tables 1-2). The deviation of 
coefficient a from 1 is caused by a small yet flnite inertia 



of the particles and by the dependence of coefficient a 
on the mean temperature gradient. The exact value of 
parameter a for inertial particles cannot be found within 
the framework of the theory of turbulent thermal diffu- 
sion (Elperin et al. 1996; 1997; 1998; 2000b; 2001) for 
the conditions of our experiments (i.e., for strong mean 
temperature gradients). However, in all the present ex- 
periments performed for different ranges of parameters 
and different directions of a mean temperature gradient, 
the coefficient of turbulent thermal diffusion a was more 
than 1, which does agree with the theory. In particular, 
the theory predicts a = 1 for noninertial particles, and 
a > 1 for inertial particles. As can be seen in Tables 
1-2, parameter a slightly decreases when the Reynolds 
number increases. 

There are other factors that can affect the spatial dis- 
tribution of particles. The contribution of the mean flow 
to the spatial distribution of particles is negligibly small. 
Indeed, the normalized distribution of the scattered light 
intensity measured in the different vertical strips in the 
regions where the mean flow velocity and the coefficient of 
turbulent diffusion vary strongly are practically identical 
(the difference being only about 1 %). Due to the effect 
of turbulent thermal diffusion, particles are redistributed 
in the vertical direction in the chamber. In turbulent 
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TABLE I: Parameters of turbulence and the turbulent ther- 
mal diffusion coefficient for a turbulent flow with an upward 
mean temperature gradient. Here Re = L (u^) jv, L is the 
integral scale of turbulence, and v is the kinematic viscosity. 
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TABLE II: Parameters of turbulence and the turbulent ther- 
mal diffusion coefficient for a turbulent ffow with a downward 
mean temperature gradient. 



the particle collision rate using the Saffman- Turner for- 
mula (SafFman and Turner 1956) confirmed this finding. 
Consequently, we did not observe any coalescence of par- 
ticles. The effect of the gravitational settling of small 
particles (0.5 — 1 /im) is negligibly small (the terminal 
fall velocity of these particles being less than 0.01 cm/s). 
Thus, we may conclude that the two competitive mech- 
anisms of particle transport, i.e., mixing by turbulent 
diffusion and accumulation of particles due to turbulent 
thermal diffusion, exist simultaneously and there is a very 
small effect of gravitational settling of the particles. 

In regard to the accuracy of the performed measure- 
ments, we found that uncertainties in optics adjustment 
and errors in measuring a CCD image background value 
are considerably less than the observed effect of a 10% 
change of normalized intensity of the scattered light in 
the test section in the presence of an imposed mean 
temperature gradient. In the range of the tracer con- 
centrations used in the experiments the particle-air sus- 
pension can be considered as an optically thin medium, 
from which we can infer that the intensity of the scat- 
tered light is proportional to the particle number den- 
sity. The total error in our measurements was deter- 
mined by {5N/N + 6f/f)/y/Q « 0.3% (see Eq. (0)), 
where 6T — Q.IK is the accuracy of the temperature 
measurements; 6N/N = 0.8% is the accuracy of the 
mean number density measurements; and Q = 8 is the 
number of experiments performed for each direction of 
the mean temperature gradient and for each value of the 
frequency of the grid oscillations. The total variation of 
the normalized particle number density due to the ef- 
fect of turbulent thermal diffusion is more than 10% (see 
Figs. 4-6). The relative error of the measurements is less 
than 3% of the total variation of the particle number 
density. Thus, the accuracy of these measurements is 
considerably higher than the magnitude of the observed 
effect. As such, our experiments confirm the existence 
of an effect of turbulent thermal diffusion, as predicted 
theoretically by Elpcrin ct al. (1996; 1997). 



flows with an upward mean temperature gradient parti- 
cles accumulated in the lower part of the chamber, and in 
flows with a downward mean temperature gradient par- 
ticles accumulated in the vicinity of the top wall of the 
chamber, i.e., in regions with a minimum mean temper- 
ature. The spatial-temporal evolution of the normalized 
number density of particles Nr is governed by the conser- 
vation law of the total number of particles (see Eqs. (0) 
and (PJ). Some fraction of particles sticks to the walls of 
the chamber, and the total number of particles without 
feeding fresh smoke slowly decreases. The characteristic 
time of this decrease is about 15 minutes. However, the 
spatial distribution of the normalized number density of 
particles does not change over time. 

The number density of particles in our experiments was 
of the order of 10^" particles per cubic meter. Therefore, 
the distance between particles is about 1 mm, and their 
collision rate is negligibly small. Indeed, calculation of 



VI. CONCLUSIONS 

A new phenomenon of turbulent thermal diffusion has 
been detected experimentally in turbulence generated by 
oscillating grids with an imposed vertical mean tempera- 
ture gradient in air flow. This phenomenon implies that 
there exists an additional mean flux of particles in the 
direction opposite to that of the mean temperature gra- 
dient, that results in the formation of large-scale inhomo- 
geneities in the spatial distribution of particles. The par- 
ticles accumulated in the vicinity of the minimum mean 
fluid temperature. In the experiments in two directions 
of the vertical mean temperature gradient, it was found 
that in a flow with a downward mean temperature gra- 
dient particles accumulate in the vicinity of the top wall 
of the chamber. In a flow with an upward mean temper- 
ature gradient particles accumulate in the vicinity of the 



8 



bottom wall of the chamber. Turbulent thermal diffusion 
can explain the large-scale aerosol layers that form inside 
atmospheric temperature inversions. 
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